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Dedication. 

During our studies in field theory we many times returned to the set of methods 
and ideas, which was developed by Misha Marinov. Here also belongs the path 
integral method in its different forms, with spin included in the quantum- 
mechanical world-line version, which is described below. The exccllenLreview 
written by Misha on the subject which appeared in Physics Reportslj was a 
most-read source on the subject at that time. It is a great pleasure and honour 
for us to dedicate the paper to his memory. 

1 Introduction. 

The present stage of development of field theory in general and of QCD in par- 
ticular requires the exploiting of nonperturbativc methods in addition to sum- 
ming up perturbative series. This calls for specific methods where dependence 
on vacuum fields can be made simple and explicit. A good example is pro- 
vided by the so-called Fock-Feynman-Schwinger Representation (FSR) based 
on the Fock-Schwinger proper time and Feynman path integral formalismcrB 
For QED asymptotic estimates the FSR was exploited in Ref. 4. Later on this 
formalism was rederived in Ref. 5 for scalar quarks in QCD and used in the 
framework of the stochastic background method inQ 

More recently some modification of the FSR was suggested in Refs. 7,8. 
The one-loop perturbative amplitudes are especially convenient for FSR. These 
amplitudes were extensively studied in Refs. 9-11. MeanwhikLthc first exten- 
sion of FSR to nonzero temperature field theory was done inE30 This formed 
the basis of a systematic study of the role nfjipnperturbative (NP) configura- 
tions in the temperature phase transitions ! 1 

One of the most important advantages of the FSR is that it allows to 
reduce physical amplitudes to weighted integrals of averaged Wilson loops. 
Thus the fields (both perturbative and NP) enter only through—Wilson loops. 
For the latter case one can apply the cluster expansion methodfj which allows 
to sum up a series of approximations directly in the exponent. As a result 
one can avoicLthe summation of Feynman diagrams to get the asymptotics of 
form factorstS The role of FSR in the treatment of NP effects is more crucial. 
In this-case one can develop a powerful method,of background perturbation 
theorjt-3 treating the NP fields as a backgroundta 

In the present paper the main focus will be on perturbative QCD, with 
the aim of establishing the correspondence between the standard perturbative 
expansion and FSR based expansion, stressing the point that FSR allows to 
make exponentiation in a very simple way. Finally also the relation of FSR to 
the popular Bern-Kosower-Strassler method is discussed. 
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2 General form of FSR in QCD. 

Let us consider a scalar particle (e.g. Higgs boson) interacting with the non- 
abelian vector potential, where the Euclidean Lagrangian is given by 

L v = \\DM 2 + \™ 2 W\ 2 = \\{d, - igAM 2 + \™ 2 W\\ (1) 

Using the Fock-Schwinger proper time representation the two-point Green's 
function of ip can be written in the quenched approximation as 

/>oo 

G{x,y) = {m 2 -Dl)^ = (x\P / dse^-^y). (2) 

Jo 

To obtain the FSR for G a second step is needed. As in Ref. 2 the matrix 
element in Eq. (^) can be rewritten in the form of a path integral, treating 
s as the ordering parameter. Note the difference of the integral (|2|) from the 
case of the Abelian QED treated in Rcfs. 2,4,9: in our case is the matrix 
operator A^{x) — A^(x)T a . It does not commute for different x. Hence the 
ordering operator P in Eq. (^) . The precise meaning of P becomes more clear 
in the final form of a path integral 

G(x,y) = J ds(Dz) xy e~ K Pexp(igJ A^dz^J , (3) 

where K = m 2 s + | J Q S c?r(^ i ) 2 . In Eq. (||) the functional integral can be 
written as 

(Dz) xy * um n /_^<me: =i *(")-(-v)) (4) 

V > xy jV^oo LLJ ( 47r£ )2 J (2^-)4 V ) 

with Ne = s. The last integral in Eq. (||) ensures that the path z M (r),0 < 
r < s, starts at £ M (0) = y p and ends at z ll (s) = x^. The form of Eq. (Q) is the 
same as in the case of QED except for the ordering operator P which provides a 
precise meaning to the integral of the noncommuting matrices A^ (z\), A^ (Z2) 
etc. In the case of QCD the forms and (f|) were introduced in Refs. 5,6. 

The FSR, corresponding to a description in terms of particle dynamics is 
equivalent toiicld theory, when all the vacuum polarisation contributions are 
also included B i.e. 



! N 

^ AH 

N=0 i=l 



7j IT J~J~ J( Dz i)xxexp(-K) Pexp I igfA^dz 
Dipexp I — d 4 xL v (x) I . (5) 
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Both sides are equal to vacuum-vacuum transition amplitude in the presence 
of the external nonabelian vector field and hence to each other. For practical 
calculations proper regularization of the above equation has to be done. 
The field A p in Eq. (@) 

can be considered as a classical external field or as 
a quantum one. In the latter case the Green's functions (A.. A) induce nonlocal 
current-current interaction terms in the l.h.s. of Eq. (|^). Such terms can also 
be generated by the presence of a </?-field potential, V(|y|) in the r.h.s. of Eq. 

(§. 

The advantage of FSR in this case follows from the very clear space-time 
picture of the corresponding dynamics in terms of particle trajectories. This 
is especially important if the currents can be treated as classical or static (for 
example, in the heavy quark case). The mentioned remark on usefulness of 
the FSR (Q) becomes clear when one considers the physical amplitude, e.g. 
the Green's function of the white state tr(ip + (x)tp(x)) or its nonlocal version 
tr[ip + (x)§(x, y)f{y)], where 0>(x,y) - to be widely used in what follows - is 
the parallel transporter along some arbitrary contour C(x, y) 



$(x,y) = Pexp (igj A^z)dz h 



(6) 



One has by standard rules 

G v (x,y) = (tr [<p + (x)p(x)] tr [ip + (y)cp(y)] ) A 

dsi / ds 2 {Dz) xy {Dz') xy e~ K ~ K ' (W) A + ... (7) 
o Jo 

where dots stand for the disconnected part, (G ip (x, x)G ip (y, y))A- We have used 
the fact that the propagator for the charge-conjugated field tp + is proportional 
to &(x,y) = <&(?/, x). Therefore the ordering P must be inverted, &(x,y) — 
Pexp(ig J Apfydzp). Thus all dependence on A^ in G v is reduced to the 
Wilson loop average 



(W) A = (trPcexpigj A^di 



(8) 

A 



Here Pc is the ordering around the closed loop C passing through the points 
x and y, the loop being made of the paths M (r)a, z' (r') and to be integrated 
over. 

The FSR can also be used to describe the quark and gluon propagation. 
Similar to the QED case, the fermion (quark) Green's function in the presence 
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of an Euclidean external gluonic field can be written as 

G q (x,y) = (ip(x)ij;(y)) q = (x\(m q + D)" 1 ^) 
= (x\(m q - D)(m q - D 2 )- X \y) 



(m q -D)l ds(Dz) xy e K <S> a {x,y), 



where § a is the same as was introduced in Ref. 2 

$<j(x, y) = -Pa exp (^9 J A^dz^j P F exp J drcr^F^ 



(9) 



(10) 



and cr M „ = j-^^u - 7„7 M ), while K and (Dz) xy are defined in Eqs. (g) and 
(|]). Note that operators Pa, Pf in Eq. (10) preserve the proper ordering of 
matrices and cr^F^ respectively. Explicit examples are considered below. 

Finally we turn to the case of FSR for the valence gluon propagating in 
the background nonabelian field. Here we only quote_the result for the gluon 
Green's function in the background Feynman gaugealia We have 



G^(x,y) = (x\(D 2 x S^ - 2tgF^)- 1 \ 



(11) 



Proceeding in the same way as for quarks, one obtains the FSR for the gluon 
Green's function 



G^(x,y) 



where we have defined 



Kn = - 



dr 



ds{Dz) xy e K °^^{x,y), 



dr. 



(12) 



^^u(x,y) 



P A exp yig j A x dz x j Pf cxp \^2g J dTF ap (z(T)) 



(13) 



[At/ 



Now in the same way as is done above for scalars in Eq. ([?]) , one may consider 
a Green's function, corresponding to the physical transition amplitude from a 
white state of q±, q~2 to another white state consisting of q^, q^. It is given by 



G^g{x,y) = (G q (x,y)rG q (x,y)r - G q (x,x)TGg(y,y)T) A) 



(14) 



where T describes the interaction between the q, q pair in the meson. The first 
term on the r.h.s. of Eq. (^) can be reduced to the same form as in Eq. 
((?]) but with the Wilson loop containing ordered insertions of the operators 
a^F^ (cf. Eq. ©). 
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3 Perturbation theory in the framework of FSR. Identities and 
partial summation. 

In this section we discuss in detail how the usual results of perturbation theory 
follow from FSR. It is useful to establish such a general connection between 
the perturbation series (Feynman diagram technique) and FSR. At the same 
time the FSR presents a unique possibility to sum up Feynman diagrams in 
a very simple way, where the final result of the summation is written in an 
exponentiated wayt3ll3 This method will be discussed in the next section. 

Consider the FSR for the quark Green's function. According to (^), the 
2-nd order of perturbative expansion of Eq. ( |ITj| ) can be written as 

G g (x,y) = (m q - D) I ds dr x / dr 2 e~ K (Dz) xu d A u{Dz) uv d i v{Dz) vy 
Jo Jo Jo 

x {igA^ujii^ + g<T^F^(u)) {igA v (v)i) v + gcr\ a F\a(v)) , (15) 
where we have used the identities 

(Dz) xy = (Dz) xu ( Tl) d i u(T 1 )(Dz) uiTl)v ( T2 )d i v(T2)(Dz) viT2)y , (16) 

pOO pS pTl pOO pOO pOO 

ds dn / dr 2 f(s, n, r 2 ) = ds dn dr 2 f(s + n + t 2 ,t 1 + t 2 ,t 2 ). (17) 
Jo Jo Jo Jo Jo Jo 

One can also expand only in the color magnetic moment interaction (<jF). This 
is useful when the spin-dependent interaction can be treated perturbatively, as 
it is in most cases for mesons and baryons (exclusions are Goldstone bosons 
and nucleons, where spin interaction is very important and interconnected with 
chiral dynamics). In this case one obtains to the second order in (uF) 

/•oo p oo /*oo 

G q 2) (x,y) = i(m q -b) ds drJ dr 2 e~ m > +Tl+T2 ^ K ^ Kl ~ K2 
Jo Jo Jo 

(Dz) xu Hx, «) 3 (^(M))(i 4 !i(Oz) u „$(u, v)g(aF(v))d i v(Dz) vy . (18) 
In another way it can be written as 

Gf\x, y) = i(m q - D)(m 2 q - Dj)" 1 ^ g(aF(u))(m 2 - tiftztfv 

xg(aF{v))(m 2 q ~Dl)^. (19) 

Here (m 2 — -D^)" 1 is the Green's function of a scalar quark in the external 
gluonic field A^. This type of expansion is useful also for the study of small- a; 
behavior of static potential, since the correlator (aF(u)aF(v)) plays an im- 
portant role there. 



7 



However, in establishing the general connection between perturbative ex- 
pansion for Green's functions in FSR and expansions of exponential $> a in Eq. 
( |l0| ) , one encounters a technical difficulty since the coupling constant g enters 
in three different ways in FSR: 

f . in the factor (m q — D) in front of the integral in Eq. (^) 

2. in the parallel transporter (the first exponential in Eq. ([To|)) 

3. in the exponential of g(aF). 

Therefore it is useful to compare the two expansions in the operator form, 
the standard one 

(m + by 1 = (m + d- igAy 1 = (m + By 1 + (m + By 1 igA{m + + 
+ (m + B)-HgA(m + By 1 igA{m + By 1 + ... (20) 

and the FSR 

OO 

(m + by 1 = (m - b)(m 2 - d 2 )- 1 ^2(6(m 2 - d 2 y 1 y , (21) 

n—o 

where we have introduced 

8 = -ig{Ad + dA)-g 2 A 2 = b 2 -d 2 . (22) 
To see how the expansion © works, using D — d — igA Eq. j2l| ) becomes 

oo 

{m + by 1 = [(m + dy 1 + igA(m 2 -d 2 y 1 ]J2[5(m 2 -d 2 y 1 ] n 

n=0 

Separating out the first term we may rewrite this as 

oo 

(m + by 1 = (m + dy 1 + (m + dy 1 igA(m- D){m 2 -d 2 y 1 ^2[S{m 2 -d 2 )- 1 }" 

n=0 

(23) 

The last three factors in Eq. ( |23| ) are the same as occurring in Eq. (pi|). 
As a consequence the formal iteration of the resulting equation for the Greens' 
function reproduces the same series as in Eq. (^0|) , showing the equivalence of 
the two expansions. 

It is important to note that each term in the expansion in powers of 5, 
after transforming the operator form of Eq. (^) into the integral form of 
FSR, becomes an expansion of the exponential <j><j in Eq. (^) in powers of g. 
The second order term of this expansion was written down before in Eq. (£L5|) . 
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It is our purpose now to establish the connection between the expansion 
@, @) and the expansion of the exponential $ CT in Eq. ( [To| ) in the quark 
propagator (|^) . One can start with term linear in A and write (for the Abelian 
case see Appendix B of Ref. 19) 

= ig J Gf\x,z(n))d i z^-A li (r 1 )G^(z(T 1 ),y), (24) 
where the notation is clear from the general representation of G q , given by Eq. 

(D 



q en 



E 



, 4e 

n— 1 



AT 

n ( ,.. - 

n— 1 x ' 

with f(n) = z(ra) - z(n - 1), ^(n) = ±L4 M (z(ra)) + A M (z(n - 1))] and 



^ r {A,^=Pexp{igJ2^{n)^{n)+gJ2^F^(z(n))e}. (26) 

n— 1 n— 1 

Representing £(ra) in Eq. (pij ) as + ^(R)), where £ M (£) refers to the 

integral over £ M in Gq ^ to the left of £ M in Eq. ( pi]) and to the integral 

in Gq standing to the right of we obtain 



- -i^-^ = 2ip^ B . (27) 
(47re) z dp M 



Thus Eq. (|24|) in momentum space becomes 

= -gGf\q)( q \p^ + A^p^Gfiq') (28) 
In a similar way the second order term from the coinciding arguments yields 

G (2) (coinc) = -g 2 [ G g a) (x, z)A 2 (z)d 4 zG Q 0) (z,y). (29) 



Finally, the first order expansion of the term a^F^ in Eq. ([Toj ) yields the 
remaining missing component of the combination S, Eq. (|22]), which can be 
rewritten as 

S = -ig(A^ + d^) - g 2 Al + ga^F^. (30) 
Hence the second term in the expansion ( pl| ) 

{m + D)- 1 = (m-b)(m 2 ~d 2 )- 1 + (m-b){m 2 -d 2 )- 1 S{m 2 -d 2 )- 1 + ... (31) 
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is exactly reproduced by the expansion of the FSR (Q), where in the first 
exponential & a in Eq. one keeps terms of the first and second order, 
0{gA ll ) and 0((gA^) 2 ), while in the second exponential one keeps only the 
first order term 0{ga^ lv F tlv ). It is easy to see that this rule can be generalized 
to higher orders of the expansion in 5 in Eq. ( pl| ) as well. 

4 Summing up leading perturbative contributions in FSR. Sudakov 
asymptotics. 

Consider now the n-point Green's function with external momenta pi at the 
i-th vertex 

G(pi, ... Pn ) =< Ji(pi)... J„(p„) >, Ji(x) = ^(x^T^(x) (32) 
Insertion of Eq. (§) into E Q- © for thc 

one-fermion loop yields 

n poo 

G(p 1 ,...p n ) =<trT\r i (m i - A) / dst 

i=i J ° 

x(Dz( , >) l( „ )I (,-i ) e- x '$i ,) e"' ( ' ,I " l (ti (i) > A (33) 

We shall disregard in what follows the factors Yi(m,i — Di) since we shall 
be interested only in the exponentiated contributions. Performing the dx^ 
integrals, one obtains 

n 

i=i 

where 

«. - / ^fi*. n < r,K - > 

(34) 

with 

n 

<w a >=< JJ$« > A , 

where is the momemtum of the fermion loop going from the i-th to the 
(i+l)-th vertex. The integral d 4 q denotes the integral over one of ; all others 
being expressed through it and all pi. 

We note that < W a > is a gauge invariant quantity summing all the 
perturbative exchanges inside the fixed Wilson contour, defined by the set 
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{£(&)}• In addition to the usual Wilson (charge) vertices, there are also 
magnetic moment vertices aF, hence the notation < W a >. 

We concentrate now on the contribution of the field in (Q), yielding 
the dominant contribution in the asymptotics (the reader is referred for the 
discussion of the aF term to the Appendix of Ref. 16. ._. 

The crucial step for what follows is the use of the cluster expansion methodjLa 
which yields for < W a >— >< W > 



(ig) r 

rl 

r=l fe 



< W >= exp{5^^-5^f w (fci)^(fej)...fMr(*r) < A^ 1 (z kl )...A t , r (z kr ) >} 



n (35) 

Here double brackets denote cumulantsEj The lowest order contribution 
(in the exponent) can be expressed through the photon (gluon) propagator. 
In the Feynman gauge it is (the gauge is irrelevant since < W > is gauge 
invariant) 

<M*)M*) >= t^{z-% - (36) 

Here C%{f) is the quadratic Casimir operator for the fundamental representa- 
tion, 1 is the unit color matrix. For QED one should replace C2I — + 1. 

Eq. ( |35|) represents the perturbative sum in the exponent, which by itself is 
an important advantage, since each term of perturbative expansion is already 
exponentiated. This property of exponentiation is well known for the Wilson 
loop, without using the powerful cluster expansion technique. In particular, 
the static QQ potential exponentiates and it can be defined through the Wilson 
loop as follows 

V(R) = -\boa o ±lnW(R,T). (37) 

In our case of possibly light quarks with masses one should carry out 
an additional integration over quark trajectories (i.e. over d£ (i) (fc)). In this 
section we consider the asymptotics of the amplitude ( |34| ) and therefore use the 
stationary point analysis for trajectories {C (^)}- Instead in the next section 
another method will be exploited, which was used for partial summation of 
perturbative diagrams in Refs. 9-11. In that section we shall also show how to 
modify this method to include also nonperturbative contributions. 



In what follows we confine ourselves to the lowest contribution (36) in Eq. 
( |35| ) and show that it yields the double logarithmic asymptotics. First of all 
one can persuade oneself that the approximation (|3^) yields in Eq. ( |34| ) all 
diagrams with exchanges of photon/gluon lines between fermion lines, all or- 
derings of lines included. For QCD this means the following: all orderings, i.e. 
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all intersection of gluon lines in space-time arc included, except that the color 
ordering of operators t a is kept fixed. Since the commutator of any two genera- 



tors t a is subleading at large N c , it means that Eq. (36) sums up all exchanges 



including intersections of gluonic lines in the leading N c approximation 



Our next point is the integration over d£(k) in Eq. (|34 ) which is Gaussian 
in the main term Ki, defining the measure of integration. Therefore we can 
carry it out by expanding the exponent in Eq. ( |34| ) around the stable fixed 
point £, which is obtained by differentiating the exponent in ( |34| ) with respect 
to £ (l) (fc). One has 



4tt 2 ^ (zW(fc) - zJ(k')) 2 

4tt 2 ^ ^ (z^(m) - z0)(n')) 4 

3,*' m>fc v v ; v 

Here e.g. z^(k) = E^iE^i^'Mi where we have chosen as the origin 
the coordinate of the first vertex, and all other coordinates are calculated 

using the connection a;W — a;^ 1 ) = X^fcLi^^W w hh the cyclic condition 
j.Cn+1) _ 2.(1) _ ' 

One can solve the system of equations (|38|) iteratively expanding in powers 
of g 2 . The first two terms are given in where one should replace £^ inside 
the curly brackets by 2iSiq^ . If one represents the exponential appearing in 
Eq. ( |3"4| ) after insertion of Eq. ( |35| ) as exp (—/(£, g)), then one can write 

' i,k 

9 2 C2(f) \- f (fc)^(fcQ , n , 4^ ,, Q x 

The Gaussian integration in (|34|) finally yields 

A^n^ e " /( ^^ trinv ' (40) 



where the matrix ip is 



ii 1 <9 2 



(41) 
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The most important for what follows is the term /(£,<?) which can be 
written as (at this point we reestablish Minkowskian metric) 

where we have defined Tj = fcj, and 

i-i 

Aij =J2 s kQ (k \i<j- (43) 

The integral in the last term on the r.h.s. of Eq. ( fiHf ) can be written as 
SiSjiq^q^Iijis, q), where 

UA^)=f f -, (i) ^ a „ ■ (44) 

The diagonal terms, 7y, with i — j do not contribute to the asymptotics and 
contain only selfenergy divergencies, which are of no interest to us in what 
follows. Therefore we shall consider only the nondiagonal terms with i =/= j. 

Let us first study the term with i = j — 1 ("the dressed vertex contribu- 
tion") and Ai,i + i = stq^ . Then Eq. (Q) is reduced to the form which will be 
studied below 



ha+i(s,q) = 



dad(3 



o Jo (asiqW + f3 Sj q(i)) 2 
dadp 



/ J (a 2 a 2 + f3 2 b 2 + 2af3(ab)) ^ 

with a = Siq^\b = s,jq^\j = i + 1. As it stands the integral ( fl5| ) diverges 
at small a,/3 (or at small r^Tj in Eq. ([42|)). The origin of this divergence 
becomes physically clear, when one expresses the distance zW from the vertex 
position, (we go over to the Minkowskian space-time) 

z^=2q^T l7 =2qW Tj . (46) 

The quasiclassical motion ( p6[ ) cannot be true for small r^, when quantum 
fluctuations wash out the straight-line trajectories. The lower limit r m i„ can 
be obtained from the quantum uncertainty principle 

AzAq ~ - z W )(q W - q U) ) ~ 1 (47) 
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Furthermore, we shall be interested in the kinematical region, where 

| ? « g Ci)| » {q^) 2 ,{q^f. (48) 
The value of T m i n then is found from ([l7]) to be 

' (49) 



Using Eq. ([19]) one can easily calculate the integral (j45|), since the term 2a/3(ab) 
in the denominator of the integrand in Eq. (|45|) always dominates. The result 
is 

T * = o MO (l+m M% ( V* +1) )^)M2(g ( V t+1 Wi)- (50) 

The integration of the general term with j ^= i — 1, i + 1 can be done using 
the expressions for the Spence functions. However in the general case the lower 
limit T m i n is inessential and the double logarithmic situation does not occur 
unless there is a large ratio, | \ \ ^> 1. 

We start with the open triangle, corresponding to the Sudakov vertex 
function asymptotics, i.e. when the fermion loop is not closed. We have 

G 3 = (-iq + m)- x T{q, q'){-iq' + m)" 1 . (51) 

In this case there is no integration over d 4 q in Eq. ( p4| ) and only one integral 
1x2 is present in Eqs. and ( pp[ ) (we disregard as before the selfinteracting 
pieces In , which do not contribute to the asymptotics) . 



Inserting Eq. (|50| ) into Eq. (|4J) and integration over ds\ds2 in Eq. ( |40[ ) 



yields the leading contribution, where in Eq. (|5C| ) the arguments Sj are taken 
to be Si = ^j. For the case of QED (C2 = 1) we find 



rfe ? ')~exp(-fln^ln«), (52) 



which coincides with the known Sudakov asymptotics. 

We turn now to the case of QCD, where the basic triangle diagram is closed 
due to color gauge invariance and try to find out whether the kinematical region 
( [48| ) plays an important role in the integral over d 4 q in Eq. ([40|). 

In the general case, when all qi are unconstrained and expressed through 
three external momenta p\ , P2 , P3 and one integration variable, the region ([i"8|), 
yielding double logarithmic asymptotics (DLA) ( |50|) , is suppressed due to large 
values of /(£, q) in the exponent. As a result the integral over dq does not lead 
to the DLA form for G3. 
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The situation changes however, if one considers instead of G3 the form 
factor, i.e. when the pole terms are factored out from the vertices 2 and 3 and 
the vertex functions appear there. To simplify matter, one can consider for 
the form factor the same representation (fiof). Under the integral one has the 
vertex functions ipi(ki),i — 1,3, where 



(pd); 



2 



WW ,_ P <3,(£!!W. (53) 

The definition ( |53| ) yields in the cm. system of particle 1 or 3 the familiar 
relative momentum of two emitted fermions. The presence of ipi imposes a 
restriction on the momenta g^, namely 

k\,kl<K 2 (54) 

where n 2 is some hadronic scale. 

Let us define in the Breit system the momenta 

pW, P ( 2 ), P ( 3 ) = ( P o,-|),(o,g),bo,f), 



where Q 2 ^> n 2 and Pq — M 2 + One can then easily see, that Eq. (54) 
constrains the region of integration over dq = d 4 qi to the region \q\ ~ k, |cfo — 



Po\ ~ K, and the conditions 



48) are satisfied. Hence in this case one recovers 



the Sudakov asymptotics (|52), where 



\qq'\^\ q (1) q {2) \~^-, q 2 ~<l' 2 ~K 2 , a -> a s C 2 . (55) 



5 Explicit path integration in FSR. Connection to the Bern-Kosower- 
Strassler method. 

The path integral (0) can be performed in a direct way. After the expansion 
of Eq. (^) in powers of and F^ v and exploiting Fourier transform for the 
latter one obtains Gaussian integrals which can be easily done. 

This procedure is similar to the method introduced in Refs. 9-11 for one- 
loop diagrams and the effective action in QED and QCD. To this end we 
extend the method of Refs. 9-11 to the case of meson and glueball Green's 
functions. We will show, that the method used in Refs. 9-11 can be simplified 
and generalized. 
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Let us start with three typical expressions in the order of increasing com- 
plexity. We may consider the one-loop effective action 

T{A} = ^ ^( s ) e " m2s -U S k >\Dz) xx trW{A tl \ (56) 
Jo s 

where £(s) is a regularization. Other expressions, which may be studied, are 
the heavy-light quark-antiquark Green's function 

G HL (Q) = rctae-'-iE^^-^E^^) n ^(W(A,)). 
Jo ^ (4^e) 2 

(57) 

or the Green's function for the arbitrary mass quark-antiquark system 



ni — 1 v ' n2 — 1 

Here we have introduced JViei = si,N 2 £2 — s 2 - 

Let us consider the expression (|57|), since Eq. ( |56|) can be managed easily 
in the same way. To proceed one expands in (Eq. |57D (W{A^)) in powers of 
A^ and uses the Fourier transform A^k) for the latter. To the second order 
in g one obtains 

Ghl(Q) - Gf L {Q) = -g 2 f W^-^M-iE^ 

Jo 



x TT 7T^-k I z»{ti)dh [ z v {h)dt2{A„{k 1 )A v {k 2 )) 



(2tt) 8 

Defining similarly to Ref. 11 

2 



y^ik^ Z {ti)=i f J^T) Z ^T)dT, J fl (T)=k^6(T-t 1 )+k^6(T-t 2 ) ) (60) 
<=1 ^ 



1G 



one arrives at the integral 



z—l n n 



fc=l 



(61) 

where we have used that z^(t) = ^2k=i£p-(k) ano - interchanged the order of 
summation in J J^(T)2 M (r)dT = J A1 (n)z At (n)e. 

The Gaussian integration in Eq. ( |6l| ) can be performed trivially and one 
obtains 



N N N 



N N 



i(q, j) = ex P {-e 3 E E - Q 2s + £22 ^ E E J »> 



fc— 1 n— /c ? 



k—1 n—k 



= exp{- / dr(/ J M (r')dr') 2 - Q 2 s + 2Q p / dr / J M (r')dr'}. (62) 
Substituting ( |60| ) into (|62|) one finally gets 
7(Q,J) = exp{-Q 2 s + 2Q At (fcWt 1 + fc( 2 )t 2 )- £ fc<?A#Wti,*j)}, (63) 

t ,.7=1,2 

where 

g B (U,t d ) = Min(t u tj). (64) 

In the derivation of Eq. (|6^) one could apply the string vertex technique of 
Refs. 9-11. A corresponding method for the Dirac propagator was used in Ref. 
20. 

So far we have neglected the factors z^{ti)z v {t2) in (|59|). To take them into 
account one could exploit the trick suggested in Ref. 11. Instead we shall use 
below another approach, based on the following relation ( we denote I^(Q, J) 
the integral (61) with z^z v taken into account). Noting that zAt) — 



-(*) 



obtain 



i„aq,j) = n 



d^jn) g M (fci) &(fc 2 ) 
L (47re) 2 e e 



x exp 



Lfc=i 



4e 



(65) 



where j M (fc) = J2n=k J^i n )- Hence one can make a simple connection 
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Computing (|S6j) one should take into account that t% < t\. In doing so one 
obtains 



i« -> -2i(Q, - fc«); if -> -2*^ - feW - fcW). (67) 
The iJL Green's function to the second order 0(g 2 ) has the form 

Ghl(Q) - Gf L {Q) = V r ifae-C^)' d4fc |^y 2) (I,(fcW)i,(fc( 2 ))) 

JO l Z7r 7 



Jo Jo 



(68) 



In some cases another form is convenient with integrals over dtidt 2 taken ex- 
plicitly, 



G HL (Q)-G%\{Q)=lg 



d 4 k^d 4 kW (i M (fcW)i„(fc( 2 >)) 
(2tt) 8 2Qfc( 2 ) - (fc( 2 )) 2 -2fc«fc( 2 ) 

1 



(Q/, ^(Q* ^ { 2 Q(Jfe(l) + jb(3))_(fe(l)+fe(2)) 



m 2 + Q 2 



m 2 + (Q-fcW-fcO 2 )) 2 



1 



1 



1 



i 2 + Q 2 m 2 + (g-fc( 1 )) 5 



(69) 

Eqs. (^8|) and ( |69| ) contain the general result for the contribution of the 
second order correlator to Gjjl(Q)- One can further specify the correlator 
{A ll {k^ 1 ')A v {k^ 2 ')) . For example, for the perturbative propagator in the Feyn- 
man gauge one has 



(A^k^)A u (k^)) = 



4ir6, 



(fc(l)_fc(2))2' 



(70) 



while for the nonperturbative correlator in the Balitsky gauge one obtains, 
using results fromEil 



(I M (fc)i,(fc')) = 



-(fe 4 -fei) 2 T„ 2 ,1 



o Jo 



dt / d£ a,j,(t)a v (t') 



d d 

xiS^Sik-Stohp) ( i-Q^i-QjJ ) exp{-[(a + f 2 )k 2 + (a + P)k' 2 
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+ 2tt'(k • k')][4a(a + P + t 12 )]- 1 } ■ [a(a + P + t 12 )]-^ 2 . (71) 

Here we have introduced above the regularizing factor exp(— a(z 2 + w 2 )) and 
used the notation 

t ~ t' 

t = = — ; Oiftit) : a 4 {t) = 1, a^t) =t,i = 1,2,3. 

11 g 11 g 

It is assumed that the limit a — ► is taken at the end of the calculations. 

Inserting Eq. (fn]) into Eqs. J6S| ) and ( |69| ) one obtains the lowest order 
nonperturbative contribution to Ghl, which can be studied for all values of Q. 
In particular, one can expand the result in powers of 1/Q 2 to compare with the 

standard OPE result, where the expansion starts with ' r(F ^ 4 (0)) . The kernel 
( |7l| ) is nonlocal and generates a series in powers of (QT g ) 2 . See Ref. 22 for 
further details. 

The equations ( p9[j69| ) obtained above are valid in the second order of ex- 
pansion of (trW(A fl )) in gA^. The method used above can easily be generalized 
to higher orders. Indeed, writing 

(trW(A^} = J2(ig) n (tr^J^ A^k^z^dtt ^ A„ 2 (fc (2) )i Ma (fc)tfta 
A, n (fcW)^ (t n )dt n e ^X*-) [J |^ } V (72) 



Tl=0 
t»-l 



II 



one can define similarly to Eq. (|6C|) for the n-th term in Eq. (72) 

n 

J^ ) (r)=J2k^S(r-t i ). (73) 



Doing the integrals over rifc=i d£(k) as before (cf. Eq. (|6l|)) one obtains the 
following result, generalizing Eq. (62), 

I< n XQ,J) = exp(-£dT Qf j M (r')dr^ - Q 2 S + 2Q M ^ dr^ J M (r')dr'| . 

. . (74) 

Inclusion of the terms i J u 1 (ii)-"i J u„(in) can be done in the same way as in Eq. 
( |66| ) and the generalization of I^ n \Q,J) to lj^...^ n (Q, J) is straightforward. 
In this way one obtains the general form of the expansion of Green's functions 
( p6| ) -(|5S|) to all orders in g. In contrast to the results of Refs. 9-11 one keeps 
not only perturbative but also nonperturbative contributions. 
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Although we have used here the FSR method to study the case of scalar 
quarks, it can also also be applied for quarks with spin starting from Eq. (p"8|). 
We have shown in the present paper how to use the FSR to treat perturbative 
QCD, establishing in particular the correspondence between the standard and 
FSR perturbation expansion. 
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